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Abstract: We find the DIS structure functions at strong coupling by calculating i?-current cor- 

j> = relators on a finite-size shock wave using AdS/CFT correspondence. We improve on the existing 

j^ ■ results in the literature by going beyond the eikonal approximation for the two lowest orders in 

CO ■ graviton exchanges. We argue that since the eikonal approximation at strong coupling resums inte- 

. : ger powers of 1/x (with x the Bjorken-x variable), the non-eikonal corrections bringing in positive 

^ . integer powers of x can not be neglected in the small-x limit, as the non-eikonal order-x correction 

O to the (n + l)st term in the eikonal series is of the same order in x as the nth eikonal term in 

that series. We demonstrate that, in qualitative agreement with the earlier DIS analysis based 

"H , on calculation of the expectation value of the Wilson loop in the shock wave background using 

%j ■ AdS/CFT, after inclusion of non-eikonal corrections DIS structure functions are described by two 

. Pr , momentum scales: Q^ ~ A^ A^^^ /x and Q2 ~ A^ A^^^, where A is the typical transverse momentum 

in the shock wave and A is the atomic number if the shock wave represents a nucleus. We discuss 

possible physical meanings of the scales Qi and Q2- 
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1. Introduction 

Over the past two decades there have been significant progress in our theoretical understanding 
of the physics of parton saturation/Color Glass Condensate (CGC) [1-27]. The Jalilian-Marian- 
lancu-McLerran-Weigert-Leonidov-Kovner ( JIMWLK) [12-19] and Balitsky-Kovchegov (BK) [20- 
24] evolution equations have been constructed which unitarize the linear Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) [28,29] evolution equation for DIS on a nucleus. The phenomenological successes of 
the CGC physics in describing the data from both the Deep Inelastic Scattering (DIS) experiments 
at HERA [30-33] and from heavy ion collision experiments at RHIC [34-37] allows one to believe 
that CGC physics correctly captures some of the main features of QCD dynamics in high energy 
scattering. 

In recent years a number of research efforts have been aimed at sharpening the quantitative 
predictive power of CGC/saturation physics. Running coupling corrections to JIMWLK and BK 
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evolution equations have been calculated in [38-41] and led to a marked improvement in the agree- 
ment between CGC predictions and the experimental data [37,42]. Subleading-A'^c corrections to 
BK evolution were analyzed in [43] and were found to be very small, though some DIS observables 
were shown to be sensitive to the difference in [44]. Next-to-leading logarithmic (NLO) corrections 
to BK evolution equation have been calculated in [45,46]: the corrections were found to be in 
agreement with the NLO BFKL calculation of [47,48] and, therefore, numerically large for linear 
evolution. While it is not clear whether NLO corrections are large in the solution of the full non- 
linear NLO BK equation, since such solution is yet to be obtained, it is important to estimate the 
size of higher order corrections to the JIMWLK and BK evolution equations beyond the running 
coupling corrections found in [38-41]. The assessment of the size of higher order correction may 
happen by performing explicit higher order calculations of the BK kernel obtaining a (presumably 
numerical) solution of the BK equation at each order. 

Alternatively, to estimate the size of higher-order corrections in the extreme large-coupling limit, 
one may use the Anti-de Sitter space/conformal field theory (AdS/CFT) correspondence [49-52] to 
study DIS. Indeed AdS/CFT correspondence is a duality between A/" = 4 super Yang-Mills (SYM) 
theory and type-IIB string theory, and as such does not apply directly to QCD. Still since A/" = 4 
SYM theory is a QCD-like gauge theory, i.e., it contains gluodynamics as a part of the theory, there 
is hope that many of its qualitative features and, possibly, some quantitative ones would apply to 
QCD. Certainly, on the perturbative side, the leading-order (LO) (pure-glue) BFKL equation is 
identical in QCD and in A/" = 4 SYM, with many similarities at NLO [53] as well. 

High energy scattering in general, and DIS in particular, in the context of AdS/CFT correspon- 
dence has been studied by many groups [54-73]. In those works the pomeron intercept at large 't 
Hooft coupling \ = g^ Nc has been calculated [54,63,65,67], though some disagreement still exists 
about its precise value [74,75]. Another important quantity for elucidating higher-order corrections 
to saturation/CGC physics is the saturation scale Qg, a momentum scale below which, in the pertur- 
bative framework, the non-linear saturation effects become important [25-27]. In CGC approaches 
based on LO BK or JIMWLK evolution the saturation scale grows as an inverse power of Bjorken-x 
variable and as a power of the nuclear atomic number A for DIS on a nucleus, Q^ ~ A^^^ (l/a;)'^""***"", 
with as the strong coupling constant. In the AdS/CFT framework the saturation scale has been 
calculated in [65,66] for DIS on an infinite thermal medium with the result that Ql ~ 1/x^ at large 
A. In [65,66], following [55,56], electromagnetic current of the standard model was replaced by the 
/^-current in A/" = 4 SYM theory. The hadronic tensor of DIS was then replaced by a correlator 
of two /^-currents, which was calculated at large 't Hooft coupling using the methods of AdS/CFT 
correspondence. Generalization of the method of [65,66] to the case of DIS on a finite-size medium 
(modeling a proton or a nucleus) was done in [70,71]. The saturation scale obtained in [70,71] 
scaled as Q^ ~ A^^^/x. 

An alternative approach to DIS was suggested in [67]: in QCD it is well-known that DIS 
at small-x can be viewed as virtual photon splitting into a quark-ant i- quark pair with the pair 
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interacting with the proton/nucleus [20,27]. Since only the interaction of the quark dipole with the 
proton/nucleus is described by strong interactions, only this part of the DIS cross section can be 
strongly coupled and should be modeled using AdS/CFT. In [67] the forward scattering amplitude 
of a dipole on a nucleus has been calculated modeling the ultrarelativistic nucleus by a shock wave 
in AdS5. Expectation value of the corresponding Wilson loop in the shock wave background was 
then calculated using the AdS/CFT prescription [76]. The dipole scattering amplitude obtained 
in [67] allowed for successful descriptions of some of the HERA DIS data [68,77], albeit in a limited 
region of small photon virtuality Q^ where QCD coupling constant should be large. While the shock 
wave considered in [67] had a finite longitudinal extent, as we will show below the results of [67] 
can be easily generalized to an infinite-size shock wave, giving the saturation scale Q\ ~ 1/x^, in 
agreement with [65,66]. However, the saturation scale for the interesting and realistic case of a 
finite-size shock wave found in [67] scales as Qs ~ A^^"^ {Ijxf ~ A^/^, in disagreement with the 
results of [69-71] (though in apparent agreement with [78], where a similar method of inserting a 
fundamental string in the bulk was used, though for the purpose of jet quenching studies). 

The goal of the present paper is to attempt to reconcile the results of [67] with that of [69-71] 
and/or to elucidate the origin of the discrepancy. We will try to perform i?-current DIS calculation 
without employing the eikonal approximation used in [69-71] . Our motivation is the following. The 
eikonal series of graviton exchanges in AdS/CFT sums up powers of 1/x on the gauge theory side. If 
X is small, this is a series in powers of a large number 1/x, and, as such, is susceptible to corrections. 
Namely, order- a; non-eikonal correction to the (n -|- l)st term in the series is (l/x)'^"'"^ x x = (l/x)", 
i.e., it is of the same order as the nth term in the series. Since the coefficients in the eikonal series 
are functions of Q^, the condition of non-eikonal corrections being small translates into a bound 
on Q^ . Below we will show that the eikonal approach of [69-71] is valid only for Q^ > (Qi)^ with 
(<5i)^ ~ A^''^ jx the candidate for the saturation scale found in [69-71]. The breakdown of eikonal 
approximation is due to the presence of another scale in the problem, Qi ~ A^/^, which corresponds 
to the candidate for the saturation scale found in [67]. Our conclusion is that _R-current DIS is a 
two-scale problem and that the exact solution of the problem should determine which of the scales 
Qi and Q2 is the saturation scale in strong-coupling DIS. 

The paper is structured as follows. We start in Sec. (0) by defining all the main concepts 
and quantities used in the calculation. In Sec. (|3.1|) we construct general exact expressions for 
the hadronic tensor modeled in AdS/CFT. The expressions for two independent components of 
the hadronic tensor are given in Eqs. ( |3.21| ) and ( |3.33D below. As these expressions appear to be 
too complicated to be evaluated precisely analytically, here we first evaluate them using the eikonal 
approximation of [69-71,79] in Sec. |3.2| . In the process we find the applicability region of the eikonal 
approximation: Q^ > (Qi)^ (see Sec. |3.2.2| ). To solidify this conclusion we evaluate Eqs. ( p.21|) 



and (|3.33|) exactly order-by-order in graviton exchanges to the first non-trivial order in Sec. |3.3| 



and show explicitly when the non-eikonal corrections become comparable to the eikonal terms. We 
summarize the results of our calculations in Sec. |3]^. Finally, in Sec. ^ we conclude by outlining 
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some of the possible physical interpretations of the scales Qi ^^^ Q2- 

2. General Setup 

Our goal is to model DIS on a shock wave at strong coupling. For simplicity we will consider shock 
waves without transverse coordinate dependence in their profile. In [80], using the holographic 
renormalization [81], the geometry in AdSs dual to a relativistic nucleus in the boundary theory 
was suggested to be given by the following metric 

t2 

ds^ = -^ {-2dx+dx- + <^(x-)z^dx-^ + dxl + dz'^]. (2.1) 

Here dx^ = {dx^y + (rfx^)^ is the transverse metric and x^ = (x'^ ±x^)/\/2 where x^ is the collision 
axis. L is the radius of S5 and z is the coordinate describing the 5th dimension with the boundary of 
AdSs at 2; = 0. Eq. ( p.l|) is the solution of Einstein equations in AdSs for $(x~) being an arbitrary 
function of x' . 

According to the AdS/CFT prescription [81] the energy-momentum tensor in the boundary 
gauge theory dual to the metric ( |2.1| ) has only one non- vanishing component: 

T-{x') = ^,Hx-). (2.2) 

Thus different functions $(x^) correspond to different longitudinal profiles of the nuclear energy- 
momentum tensor. 

Here we take a shock wave made of homogeneous matter with a finite longitudinal extent [79,82] 

<l>(x-) = ^e{x-)e{a-x-). (2.3) 

a 

While this is a simple ansatz, it appears to be quite realistic for a large ultrarelativistic nucleus. 
Indeed transverse coordinate dependence is neglected in Eq. (|2.3|) , but since the relevant transverse 
distance scales for a DIS process are much shorter than the length of a typical variation of the 
nuclear profile in the transverse direction, we believe neglecting transverse dependence in Eq. (|2.3| ) 
does not affect the physics in a qualitative way and is a good first approximation for studying DIS 
in AdS/CFT, which can be easily improved upon later. 

The parameter /i is related to the large light-cone momentum of the nucleons in the nucleus 
p"*", the the atomic number A and the typical transverse momentum scale A by [82] 

/iocp+A^A^/l (2.4) 

The longitudinal width of the nucleus is Lorentz-contracted and is given approximately by [82] 

^1/3 

a ~ -^. (2.5) 

p+ 



-A- 



Following [65,66] we will model electromagnetic current in DIS by an i?-current J^{x), which 
is a conserved current corresponding to a U{1) subgroup of the SU{A) i?-symmetry of the Af = 4 
SYM theory in four dimensions. The current J^(x) can be written in terms of the scalar, spinor 
and vector fields of the A/" = 4 SYM theory. For a systematic and pedagogical definition of the 
i?-current we refer the reader to [83,84]. 

We want to calculate the retarded i?-current correlator [65,66,71,85-88] 

W{q) = —i d^x±£y±dx-dy-d{x+ -y+)e-"'-^''-y^9{x^-y^){p\[J^'{x),r{y)]\p) (2.6) 
a J 

where all integrals run from — oo to +oo. Our metric in 4-dimensions is r]^,^ = diag{—l, +1, +1, +1). 
By \p) we denote the proton or nucleus state which is modeled in AdS/CFT by the shock wave 
( p.l| ). The essential ingredient of Eq. ( |2.6| ) is the retarded Green function in the coordinate space 



W''{x,y) = teix'-y'){p\[J^^{x),.r{y)]\p). (2.7) 

Note that, as usual in DIS, current conservation and Lorentz symmetries demand that 11^'^ (g) can 
be written in the following standard form 

n^^(g) = (^v'-' - ^) ni(x, Q') +{p'-^-^ q"^ {p"" - ^ ?^) ^2{x. Q') (2.8) 

where p is the momentum of (a nucleon in) the shock wave, 

Q2 = q\ (2.9) 

and the Bjorken-x variable is 

x = —^ . (2.10) 

-2p-q ^ ' 

(Note again our ?7^j, = diag{—l, +1, +1, +1) metric convention.) The imaginary part of the corre- 
lator U^'^lq) is proportional to the DIS hadronic tensor: nevertheless, for brevity, we will refer to 
n^'^(g) itself as hadronic tensor. 

To calculate the retarded Green function at strong 't Hooft coupling using AdS/CFT corre- 
spondence one makes use of the fact that the i?-current J^ is dual to a Maxwell gauge field in the 
bulk [86-90]. The action of the Maxwell gauge field in empty AdSs space and in the space described 
by the metric (|2.1| ) is 



^Maxwell = 'J^ [ d'x ^^ FmN F^'' = -^ [ d'x ^ FmN F^'' . (2.11) 

64 TT"' L J 64 TT^ J z'^ 

Here and throughout the paper indices M, N run from to 4, while /i, v run from to 3. 



Classical sourceless Maxwell equations in the curved background read 

dM[V^9 9''''9'''F^s\ =0. 
In AdSs the classical Maxwell action can be written with the help of Eq. ( ^.121 ) as [66] 






cl 
Maxwell 



A^.' 



32 7r2 



'4~ 



d X 



- {A+d,A_ + A_d,A+ - AAAi 



(2.12) 



(2.13) 



z=0 



with i = 1,2 denoting transverse spatial dimensions and with summation assumed over repeated 
indices. In arriving at Eq. ( |2.13|) we have made use oi A^ = gauge, which we will employ from 
now on. 

In the background of the metric ( |2.1| ) with the shock wave profile ( p.3| ) Maxwell equations 
become (labeled by component) 



[+) [z d+ d^ + d, - z d"^^] A4x+,x-,z)-zdlA+{x+,x-,z) 

= ^zH(x-)e(a-x-) \3d, + zdl] A+h 
a '- -' 

[-) [zd+d^ + d,-zdl] A+{x+,x-,z) 



2d.d- 



z 



z (9^ y4_(x^,a; , z) 



Ai{x^,x , z) = z* 9{x )9{a — x )9^Aj(x"'",x ,2;), 



(2.14a) 
(2.14b) 

(2.14c) 



(z) d,[d^A+{x+,x-,z) + d+A_{x^,x~,z)\ = ~- zH{x-)9{a- x-)d,d+ A+{x+ ,x~ ,z). (2.14d) 



In arriving at Eqs. ( |2.14| ) we have assumed that the gauge field A^{x,z) is independent of the 
transverse coordinates x_ = [x^, x^). The reason for this assumption will be explained later. 
The AdS/CFT prescription for calculating this retarded Green function (|2.7|) is [86-90] 



U>^^ix,y) 



"^Maxwell 

6Alix)6AM 



(2.15) 



where A^ (x) is the value of the classical Maxwell gauge field at the boundary of AdSs. However, 
the correlator one would obtain from Eq. ( |2.15| ) in the background of the metric given by Eqs. ( |2.1| ) 
and ( ^.31 ) would contain both the vacuum component and the /i-dependent term due to DIS on a 
shock wave. Since we are interested in the latter we need to subtract the vacuum piece. We thus 
write using Eq. (|2.6| ) 



m'^iq) 



A^ 



d x±d y±dx dy d{x^ — y^) e 



+\ f,-ig-{^-y) 



dA^MdAtiy) 



['-'Maxwell lA'') "^Maxwell (0) J ' 

(2.16) 



where S^^^^^i^fi) is the classical Maxwell field action in the background of the shock wave metric 



( p.l|) , while 5'^axweii(0) is the same action in the empty AdSs background. 
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3. i?-Currents Correlator 



3.1 General Expression 

3.1.1 Transverse Components of the Hadronic Tensor 

While Eqs. ( p.l4| ) are hard to solve exactly, we will look for the solution perturbatively in fi. We 
start by concentrating on the transverse field component Ai which contributes to the transverse 
part n*-' of the hadronic tensor 11^'^. Note that the equation ( |2.14c| ) for A^ completely decouples 
from the rest of Maxwell equations ( p.l4|) : therefore we can treat Ai as an independent degree of 
freedom. We write 



l(0) 



A{x,z) = Af\x,z) + A\''{x,z)+A\''{x,z) + 



m, 



l(2)/ 



(3.1) 



iW 



where the term A- is of the order fx^. 



Start by putting /i = (no shock wave) and solving Eq. ( p.l4c|) for A\ 



i(o), 



X, Z) 



2d+d- 






l(O) 



Ar{x\x-,z) 



0. 



(3.2) 



Concentrating on the 2;-dependence of A\ (x, z) we see that the general solution of this equation 
can be written as 



Af\x+,x-.z) 



^2d+d^Ki{z ^2d+d^) Ci(a;+, x') + z ^2d+d^h{z ^2d+d^) ^2(0;+, 



X 



(3.3) 



with Ci and C2 some arbitrary functions. Demanding that our Maxwell field (and, more importantly, 
its field strength) grows slower than ~2;^as^— J-oo^we can discard the second term on the right 
of Eq. (|3.3| ) since it would give an exponential divergence at large- 2; for positive eigenvalues of the 
operator 2(9_|_9_. We therefore write 



Af\x^,x-,z) = z ^2 d+ d. Ki{z ^2 d+ d.) 4^1 



X ,x 



(3.4) 



where A^{x) is the boundary value of the field Al (x, z). 

We define the Green function of the operator on the left-hand-side of Eq. ( p.l4c| ) (while treating 
it as a differential operator in z only) by 



2d+d- + -d,-d^^ 



G{z]z';d+d-) = z'6{z-z'). 



(3.5) 



The Green function G(z; z'; d+ d-) is itself a differential operator being dependent on 9+ d-. Unfor- 
tunately Eq. (|3.5|) does not uniquely define the Green function G since we can always shift the Green 



^This condition arises in deriving Eq. ( |2.13 ), where the contribution from z — 00 can be neglected only if all field 
components grow slower than ^ z^ at large-z. 
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function by the right-hand-side of Eq. ( p.3| ) unless we specify the boundary conditions. As one can 
see from Eq. ( p.l6| ) our goal is to differentiate with respect to the boundary values of the Maxwell 
field. At the same time, if we solve Maxwell equations ( p.l4| ) order-by-order in /i the boundary 
value of the free field in Eq. ( |3.4| ) may be modified by /^-dependent corrections and may become a 
^-dependent function itself. In other words the boundary value of the field Ai{x,z) would not be 
A'l{x) from Eq. (P^), but instead would contain some explicit /x-dependent terms added to it: it 
would then be unclear how to perform the functional differentiation of Eq. (|2.16| ). To avoid this 



problem it appears easiest to follow [71] and demand that Af' {x, z) is zero at 2; = for all n > 1. 
Then the boundary value of the full field Ai{x,z) would be given by A'^{x~^,x~) from Eq. (|3.4|) , 
making functional differentiation possible. 

We therefore demand that the Green function G{z; z']d+d^) is at 2; = 0. On top of that 
we demand that G{z; z';djf-d-) does not diverge exponentially as 2 —t- 00. These two conditions 
together with Eq. ( \i.b\) fix the Green function G{z; z'\ 9+ 9_) uniquely. The Green function can be 
shown to be equal to [66] 

G{z; z- d+d-) = zz h (z< ^2d+d\ Ki (zy y/2d+dJ) (3.6) 

with 

2;>(<) = max (min) {2;, 2;'}. (3.7) 

In general the convolution of Green function G{z; z'; 9+ 9_) with other functions may generate 
inverse powers of 9+9_. Requiring causality we will understand those as denoting the following 
operations 



X 



+ 



^[. . .](a:+) = I dx'^ [. . .](x'+), ^[. . .]{x-) ^ J dx'- [. . ]{x'-). (3.8) 

—00 —00 

Let us define one more abbreviated notation: 

/■°° dz' 
G,{d+d_)fix^,x-,z)^ / —Giz;z';d+d^)fix+,x-,z') (3.9) 

Jo z 

for an arbitrary function f{x~^, x~ ,z). 

With the help of this notation we write the result of the first iteration of Eq. ( p.l4c| ) as 



A^\x+,x-,z) = -^GJd+d)z^9(x~)9(a~x-)dlAf\x+,x-,z) 
a 



-^ G,{d+ d_) z^ e{x') e{a - X-) dl z y/2d+d_Ki{z ^2d+d_) A'l{x+, x"), (3.10) 



where in the last step we have used Eq. (3^). Repeating the procedure several times a general term 
in the series of Eq. ( |3.1| ) can be written as 



(3.11) 



Af'\x^,x-,z) = \-^G,{d+d-)z^9{x-)9{a-x-)dlT z^/2^^Kliz^/2^^)A. 



where each differential operator in the square brackets acts on everything to its right, that is, say, d- 
in one of the brackets acts on all ^"-dependence in all other brackets to its right and on A'^{x~^, x~). 
The solution of Eq. ( p.l4c|) can then be written as an infinite series 



A{x+,x-,z) = J2 [-- ^-(^+ 9-) z" e{x') e{a - X-) dlT z^y2d+ d- Ki{z ^2d+ dJ) A\{x^ 

n=0 



,X 



(3.12) 

In evaluating structure functions using Eq. ( |2.16| ) it will be important to know the small-0 
behavior of Ai{x~^, x~ ,z). Expanding Eq. ( p.4| ) in powers of z yields 



Af'\x+,x-,z) = {1 + 



z^ 2 94- d. 



In 



z"" 2 d+ g, 
4 



27-1 



+ o{zUnz)^Ali 



x+,x-). (3.13) 



iW 



For Al with n > 1 the expansion is different. Expanding the Green function in Eq. (^) at small-z 
we write for n > 1 



A^\x-^,x-,z) = _^2 \/2g+g, /i I dz'Ki{z'^2d+d-)z"'e{x-)e{a-x-)di 



2 a 



-^ G-Ad+ 9_) z'^ 9{x-) 0{a - X-) d: 



-\ ^2 



n— 1 



z' ,/2d^Ki{z' ^/2d^) A^(x+, X-) + o (z^) . 

(3.14) 



We are now almost ready to evaluate W^{q). Using Eq. ( p^.l3 ) in Eq. 
asymptotics found in Eqs. ( p.l3|) and (|3.14|) we obtain 



along with the 



A A^ /" /■ 

W^(q) = '— (fx^(fy±dx~dy-d(x^ -y+)e-"'-^''-y^ d\ 

a 16 11'^ J J 

16Af\x,z) 6Ak{x,z) _ 1 6Af\x,z) ^ SAf{x,z) 



z 6A'^{x) M (y) z 6A'^{x) M (y) 



(3.15) 



z=0 



where we made use of the fact that W^(q) is an even function of q^ (see Eq. ( p.8|) ). The sum over 
k runs over A; = 1, 2. 

Just like [71] we will work in a frame with q = 0. This is the reason we have neglected transverse 
coordinate dependence of the classical Maxwell fields throughout the discussion. The argument is 



as follows: since the metric tensor in Eq. (p.l|) is independent of transverse coordinates, putting 
transverse coordinate dependence back into Maxwell equations ( |2.14|) would only generate some 
differential operators dx^ in its solution ( p.l2| ). After substituting Eq. ( p.l2| ) into Eq. ( p.lSp , those 
differential operators become dx^ acting either on 5^(x — x) or on 6'^{y — x). We can therefore 
replace d^^ — )■ —d^^ (or d^^ — )• —dy^) and integrate by parts, such that in the end we would get 



dx^ -^ -dx^ 



-^ —ic[± (or dxj^ —7- —dy^ — 7- iq±). Hence these transverse coordinate derivatives 



vanish in the q± = limit, leaving only two transverse delta- functions S'^{x — x) 6'^{y — x), which 
eliminate x± and y± integrals. The remaining x± integral is strictly-speaking infinite, but we assume 
that the nucleus has a large but finite transverse extent and replace 



(fx± — > 5* J 



(3.16) 



where Sj_ is the transverse area of the nucleus. We assume that as long as the nucleus is large 
enough in transverse direction, much larger than the typical relevant distance scale for DIS, the 
DIS process would most of the time be insensitive to the edge effects justifying the approximation. 
This is done in complete analogy with perturbative DIS calculations [20,21]. We conclude that at 
g_L = all transverse integrals simply disappear, leaving only the factor of S± from Eq. (|3.16| ). 
Using Eqs. (|3.13|) and (p.l4|) we can find the functional derivatives needed in Eq. (|3.15|) : 

6Af\x+,x-,z) 



5A^{y+,y-) 



5^^ (5(x+ - y+) 5{x~ -y~) + o {z^ In z) 



(3.17) 



and for n > 1 



5A^^\x+,x-,z) 
5A\{y+,y-) 



X 



.^ik ^2 VE^±^ l^ I dz'K^{z'^/2d^)z'H{x-)e{a-x-)d^, 

Zj (Jj 



fl 



- G,id+ d^) z'^ e{x-) e{a -x-)d] 



n-l 



z ^2d+d-Ki{z' ^2d+d^] 



X 5{x'^ — y^) 6{x — y ) + o [z'^) 



(3.18) 



Here S^'^ is the Kronecker delta. 

Substituting Eqs. (|3.17|) and ( p.l8|) into Eq. ( p.l5|) and employing the arguments which led to 
Eq. ( p.l6|) we therefore get 



m^(g+,g-,g = 0) 



A^S^ Nl ^ 



-^ A 6'^ / dx- dy- d(x+ - y+) e''^^ (^"-s/" )+»'?- (^+-y+) J2d^_ 
n^ 2 a J 



X dz Ki{z ^2d+d.) z^ 9{x-) 9{a - x") dl ^ [-- G,{d+ d^) z^ e{x-) d{a - x~) & 



n=l 



X Z 



^2d+d^Ki{z ^2d+d^) 6{x^ - y+) 6{x- - 



n-\ 



(3.19) 
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Integrating by parts in Eq. ( 3.19|) we can replace all 9_)_ with 

(9+ — 7- —i q~ 



(3.20) 



after which the integral over x^ — y^ simply cancels 6{x^ — y^). For all factors of d- to the left of the 
sum in Eq. ( |3.19|) integration by parts also gives d^ — > —iq^. Finally, the factors of 9_ to the right 
of the sum in Eq. ( p.l9| ) are also replaced by d- — )■ —iq^, which can be seen by integrating over 
y~ in Eq. ( p.l9 ), applying d- which appear to the right of the sum, and undoing the |/~-integral. 
Finally remembering that Q^ = q^ = —2 q^ q^ > in DIS we write 



W^{q+,q-,q = 0) 



a 8 7r2 2 a 



^6''QHq-)^ I dx- I dy-e'^^^'^'-y'^ f dz z^ Ki{Q , 



n=l 



J2 \-iq'fGzi-iq~d^)z^e{x^)e{a-x~) z Ki{Q z) 5{x- - y-). 



n~l 



(3.21) 



We have purposefully did not carry out the y integration as the expression in the form shown in 
Eq. ( ^.211 ) will be easier to evaluate. 

3.1.2 Longitudinal Components of the Hadronic Tensor 



To find the structure functions Hi and 112 from Eq. (|2.8|) we need to find one of the longitudinal 
components (n^"*", 11 , or n+~) of the hadronic tensor U^'^ . Finding only one of the longitudinal 
components of U'^" is sufficient, along with Eq. ( |3.21| ), to uniquely determine Hi and 112. We will 
determine 11"''"*'. 

We start by solving Eq. ( p.l4b|) for A_: 



A^{x^,x , z) 



1 



Plugging this into Eq. ( |2.14d| ) we write 
1 1 



2d+d^-- + -d,-di 

Z'^ z 



d+d^ + -d,- dl 



fi 



A^{x^,x ,z). 



(3.22) 



d,A+(x+,x~,z) = ~^z^e(x-)e(a-x-)dld,A+(x+,x-,z). (3.23) 

a 



Similar to the transverse field components case above, we begin by solving Eqs. ( p.22|) and ( p.23|) 
in the /i = case of no shock wave. Solving Eq. ( p.23| ) for /x = while requiring that A\_' (x"*", x~, z) 
remains finite as 2; — )■ +cxd we get 



l(0) 



A^^\x^,x-,z) = Al{x-^,x-) + 



2d,d^ 



1 - z ^2d+d^Ki (z ^/2^+^^)] C{x^, 



(3.24) 
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with C(x+, x~) an arbitrary function of a;+ and x~ and A^(a;+, x^) the boundary value of the field 
A\_ (x^ , x~ , z) . To fix C{x^,x^) we plug Eq. ( p.24|) into Eq. (|3.22|) and match terms at z = 
obtaining 



C{x+,x-) = dlAt{x+,x-)-d+d^Al{x+,x~), 
where A'L{x~^, x~) is the boundary value of the field A_ {x^, x~ , z). We thus have 



(3.25) 



A'f>(x+,x-,z) 



A^^\x+,x-,z) = ^^ 



1 
~ 2 

^1 
9_ 2 

9_ 1 

9^2 L 

1 

+ 7: 



(X ,x 



l + z^/2d^Kiiz^2d^\ A^(5 
l-z^2d+d^Ki {z^2d+d\ A^(x+,x") 

1 - z ^2 9+ 9_ Ki (z ^2 d+ dS]] A\{x-^,x-) 
l + z^2d+d^ Ki (z ^2 d+ a_ 1 Atix+.x-). 



(3.26a) 



(3.26b) 



Using the /i-expansion technique developed above we can write the solution of Eq. ( |3.23|) for 

d, A. as 



d^A+{x^,x ,z) 



J2dzA^+\x+,x-,z) = J2 \-- G^{d+d_) z^ e{x-) eia - X-) dl 

L (Jj 
n=0 n=0 

X zKo{z^y2d+d^) [dlAtix+,x-)-d+d^AX{x+,x- 
where we have defined 

G^id+d^)fix^,x-,z) = 



dz 



'-GL{z;z';d^d^)f{x^,x-,z') 



with the longitudinal Green function defined by 



z^ z 



GLiz;z';d+d_) = z'5{z-z'). 



(3.27) 



(3.28) 



(3.29) 



Requiring that Gl{z; z'\ 9+ 9_) goes to zero as z — )■ and that it is finite at z -^ +oo one readily 
obtains [56, 66, 71] 



Gl{z; z'; d+ 9_) = zz'Io (^< y^2d7d^) K, [z^ y^2d^) . 



(3.30) 



The subscript L for Green function Gl stands for longitudinal components. (Here we are using the 
same notation as in [71].) In arriving at Eq. (|3.27|) we have again demanded that A\^ {x'^^x~^z = 
0) = for n > 1, such that A^{x^,x~ ,z = 0) = A^_^{x~^,x~). Eq. (|3.22|) can be used together with 



Eq. ( |3.27|) to find A_{x~^, x , z) as a series in powers of fi: one can easily show that A_ {x~^, x 



0) = for n > 1 as well. 
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To evaluate n++ we use Eqs. ( p.l3| ) and (|2.16| ) along with the fact that A^J^'{x^, x , z = 0) = 
and y4_ (x+, a;~, z = 0) = for n > 1 obtaining 



A2 JV2 

n++(g) = -^ ^ 



a IGtt^ 



d^x^ d^y± dx-dy- d{x+ - y+) e"* « ' (^"^^ / d^x 



4~ 



16Af\x,z) 6A_{x,z) 
~z 6A\{x) " ^A\{y) 



\bA^l\i,z) bA^{x,z) \bAf{x,z) bAi^{x,z) _ 1ML°^(£,z) Mf (5,z) 
^i bA\{x) " bA\{y) ~~z bA\{x) " ^A\{y) ' ~z bA\{x:) " bA\{y) 



z=0 

(3.31) 



With the help of Eqs. (p|), (p^), (pSJ) and (lO^) we write 



(5(x^ — y^) 5(a; — y ) + o (2;^ In z) 
o (z^ In 2;) 



^A\{y^,y-) 

bA\{y^,y') ~~d+ SA\{y^,y-) + ^^"^ ^ 



i(«)/ 



—2; 



/i 9_ 



(3.32a) 
(3.32b) 



dz' Ko{z' ^2d+d-) z'^ 9{x-) 9{a - x') d] 



n-l 



-^ Gf',(9+ 9_) z'^ 9{x-) 9{a - X-) d'i z' K^{z' ^2 3+ 9_) 9+ a 



X 6{x^ — y^) 6{x — y ) + o [z^) 



(3.32c) 



with Eq. ( ^32cD valid for n > 1. 

Using Eqs. ( p.32| ) in Eq. ( p.31| ) and integrating over £, x^, y±, and x"' 
was done in arriving at Eq. (|3.21|) yields 



y^ similar to how it 



n++(g+,g-,g = 0) 



A^S± Nl /i 
a 32 7r2 2 a 



a 00 



Q' /"dx" f dy^e^'^^^'^'^y'^ f dzz^KoizQ) 



n=l 



J2\-{(l~yGz{-M'd^)z^0{x')0{a-x-) zKo{zQ)6{x- -y-). (3.33) 



n-l 



Eqs. ( p.21|) and ( p.33|) give us the most general expressions for the two components of the 
hadronic tensor that we need to determine the structure functions Hi and 112. We will now evaluate 
them in the eikonal approximation. 
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3.2 Eikonal Approximation and Its Applicability Region 
3.2.1 Eikonal Hadronic Tensor 

Let us evaluate the expressions ( |3.21|) and ( |3.33|) in the eikonal approximation first. Eikonal ap- 



proximation corresponds to the mathematical limit of a — )■ 0, such that Eq. ( |2.3| ) becomes 



<l>(x ) -)■ /i(5(x 



(3.34) 



Eikonal approximation in AdS/CFT has been employed before in [61,64-66,69,70,79]. The eikonal 
approximation (|3.34|) tries to mimic the physical limit when the shock wave is infinitely boosted. 
Indeed for a real-life proton or nucleus the physical limit of infinite boost can be achieved by sending 
the large momentum of the proton/nucleus wave p"*" to infinity. While indeed in the p+ — )■ oo limit 
a — > as follows from Eq. ( |2.5| ), one also notices from Eq. (|2.4|) that strictly-speaking in this limit 
yU — 7- oo making Eq. ( p.34| ) meaningless. We will understand the eikonal limit as the case when p"*" 
is very large but is still finite, such that the delta-function approximation of Eq. p.34| ) is valid, 
though fi is very large. The eikonal approximation can be thought of as taking p^ -^ oo limit, while 
simultaneously sending A — )■ in such a way that /i would remain constant, the possibility of which 
follows from Eq. ( p.4|) . This is similar to Aichelburg and Sexl's construction of ultrarelativistic 
black hole metric [91]. Mathematically the eikonal limit is simply equivalent to taking a — )■ while 
keeping fj, fixed. 



X 



Z w 



gravitons 




O O 

o 

o 

o 

• • • o 

o 
o 
o 

^AAAAAAAAAAAAAAAAAKAAA/ 



Shock Wave 



Maxwell gauge field 



Figure 1: Diagrammatic representation of the Maxwell gauge field scattering on the shock wave in the 
bulk. The horizontal wiggly line represents the gauge field, while the vertical cork-screw lines represent 
graviton exchanges with the shock wave. The boundary of AdSs is at the top of the shock wave, which in 
turn is denoted by the shaded rectangle. 
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Following [79,82] one can identify the series in Eqs. ( 3.21 ) and ( 3.33|) with the scattering of 
the gauge field in the graviton field of a shock wave. Each power of /i in Eqs. ( p.21j ) and ( p.33|) 
corresponds to a graviton exchange with the shock wave.^ The scattering of the gauge field in the 
shock wave is shown in Fig. |l|. If a is small the shock wave has a very short extent in the x~- 
direction, such that the derivative d~ in Eq. ( |3.5| ) is very large. We thus approximate the eikonal 
Green function as (see [79] for a similar approach to shock wave scattering with the goal of modeling 
heavy ion collisions) 



G''\z-z'-d+d. 



2d+d^ 



z'6{z-z'). 



(3.35) 



One then has 



GTid+d^)fix'',x-,z) 



1 



2d^d- 



f{x^,x ,z). 



(3.36) 



We start by evaluating the transverse components of the hadronic tensor first. Replacing 
9+ — i- —2 q^ in Eq. (|3.36|) and using the result in Eq. (|3.21|) yields 



a2 c ]\r2 .. 



X 



E 



8 7r2 2 a 



i — q —- z* 9{x )9(a~x' 
2 a (7_ 



dx' 



dy e 



Now since for n > 2 and for a;~ < a 



— Q{x-')Q{a-x') 



n-l 



we obtain 





n-l 







zK^{Qz)b{x- ~y-). 



5(x -y ) = 9{y )e{a-y )e{x -y ) 



\x 



\n-2 



(n-2)! 



(3.37) 



(3.38) 



dx 



dy e 



i q+ {x -y ) 



— — 6(x ) 6(a — X 

O- 



n-l 



6{x 



a" >^ 

n! ^ 

m=0 

a"" 
(n-l)! 



(ig+ ay 



{n — l)n 



ml 



-,iq^ a 



{n + m — 1) {n + m) 



-iq a) 



\q+ a + n-l) \r{n) - V{n, -i g+ a)] | . (3.39) 



^An interesting (but irrelevant for presented calculations) question is where exactly the graviton field of the shock 



wave originates: since our shock wave (2.1) has no source in the bulk, one can think of is as having a source at z = oo 
(see [92]) with graviton exchanges with that source, or one may think of the boundary condition at z = (that we 
have a nucleus in four dimensions) as being the effective "source" for the shock wave, with the gravitons exchanged 
with the boundary, as shown in Fig. P. 
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The nth term in the series of Eq. ( p.37|) brings in a factor of fi"-/a'^~^^ on top of the factor we 
obtained in Eq. (|3.39| ). Note that one factor of 1/a comes from the prefactor in the definition of 
U'^'^ in Eq. ( |2.6| ). Therefore, the a — ?• eikonal hmit should not apply to this factor. For the purpose 
of the eikonal approximation the nth term in the series of Eq. ( 3.37]) is then of the order /i"/a"' 
times Eq. ( p.39| ). It is then clear that in the a — t- limit with fi fixed only the first term in the 
series in the second line of Eq. (|3.39|) survives, as 



a cc 

dx~ 



dy e 



ig+{x -y ) 







— e{x-)e{a-x 



n-l 



5{x- -y-) = - + o (a"+i) 



n\ 



(3.40) 



Using the eikonal approximation of Eq. (|3.40|) in Eq. ( p.37|) we obtain 



nl(g+,g-,g = 0) 



A^S, N} 



a »7r^ 



5'^Q\~ / dzz [K^iQzf 



- ~4 



1 - exp ( - /i g z 



(3.41) 



in agreement with the eikonal formulas used in [69, 70] and recently derived in [71] for DIS on a 
shock wave. Remembering that a = A^^^/p~^ we see that q~ /a = p^q~/A^^^ = Q"^ / {2 x A^^^) such 
that Eq. ( |3.41|) can be re-written as 



lll,{q\q-,q = 0) =zA'S_ 



2 o c rij V 



16 7r2 



X 



Al/3 



dzz [Ki{Qz)y 



1-exp I -/ig z 



(3.42) 



making the agreement with [71] manifest (up to a trivial factor of A^ S± which probably signifies a 
slightly different overall normalization used in [71]). 

Similar calculations for 11''"+ from Eq. ( p.33|) give the eikonal expression 



n+^(g+,g-,g = 0) = tA'S, 



Nl Q' 



Q' 



16 7r2 4(g-)2 x A^^ 



dzz [Ko{Qz)f 



1 - exp ( - /i g z 



(3.43) 



Eq. ( p.43|) agrees with the result of [71] up to the same overall normalization factor of A^ S^ as for 



n^^ 



eik' 



Note that the eikonal approximation developed in [79] also leads to the eikonal propagator (the 
factor in the square brackets) in Eqs. ( p.42| ) and (|3.43|) originally obtained in [69,70]. To see this 
note that the truncated eikonal graviton amplitude in Eq. (3.29) of [79] is proportional to 



n=0 ^ +^ 



^zH,{x^)^) t,{x^), 



(3.44) 
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where, for the purposes of comparing with Eq. ( p.42 ) we take the shock wave profile to be 

t^{x+) = ^e{x+)e{a-x+). 

02 



(3.45) 



In [79] a proton-nucleus collision was modeled by colliding a shock wave with a smaller energy 
density (a proton) with a shock wave with a larger energy density (a nucleus). The larger "nucleus" 
shock wave was chosen to move in the x~ direction in [79]: this is why its profile in Eq. (|3.45|) is a 
function of x"*" instead of x~ . 

Using Eq. ( p.45|) one can readily show that 



lim t2(x+)— t2(x+) 



^'t' 



■6{x+), 



which, when substituted into Eq. ( p.44| ) after summing the series over n yields 

1 
exp 



(3.46) 



z ll2d-\ t'i{x 



(3.47) 



Identifying c}_ in Eq. ( |3.47D with —i q^ in Eqs. ( pj.42[ ) and 



we see that the exponents in both 
equations are identical up to x^ <-)■ x~ interchange. (In [79] the series started from n = since 
n = 0-term corresponded to no additional rescatterings, but still contained graviton production, 
in contrast to the DIS case at hand, where no rescattering implies no interactions and hence no 
contribution to DIS cross section. This is why we do not subtract 1 from the exponent in Eq. ( p.47|) 
unlike the exponents in Eqs. ( p.42 ) and ( 3.43|) .) 

Finally, as /x = p^ A^ A^^^ we see that fiq~ = Q^ Is? A^f^ /(2 x). Defining a momentum scale^ 



Ql{x,A) 



AMV3 
4x 



(3.48) 



we recast Eqs. ( p.42| ) and ( |3.43| ) into the following form: 
nl(9+,r,i = 0) =iI^^S^^^J^^-^^ I dzz[K,{Qz)f [l~e^^{tQ'Ql{x,A)z^)], 



16 7r2 



X 



AV3 



n+^(g+,g-,g = 0) =iK'S^ 



m 



Q' Q' 



16 7r2 4(g-)2 x A^/^ 



(3.49a) 
dz z [KoiQ z)f [1 - exp (z Q^ Ql{x, A) z')] . 

(3.49b) 



•^As one can see from Eqs. (2.4) and ( p^ ) the expressions we use for /i and a are accurate up to a constant and 
possibly some factors of 't Hooft coupling A [82]. Such factors, which are important for phenomenology, can be easily 
included later. 



17 



In [69-71] the scale Qf was identified with the saturation scale Q^ as in the eikonal approximation 
it is the only momentum scale in the problem and because structure functions become independent 
of X for Q^ < Ql, in agreement with what one observes in perturbative approaches [21]. Our 
understanding of the physical meaning of Ql will be detailed in Sec. ^. 

3.2.2 Applicability Region of the Eikonal Approximation 

The question we would like to address now is whether the eikonal result (|3.49|) gives us the complete 
"hadronic tensor" U^'^ in the small-x limit. As we noted above, the proper physical high energy limit 
corresponds to increasing the proton/nucleon momentum p"*", and not to simply taking a — )■ limit. 
In arriving at Eqs. ( p. 491) we have made several approximations. In particular we have neglected 
higher powers of q~^ a in approximating Eq. (|3.39|) with Eq. ( |3.40|) . Since 



q a 



p+ 



-X 



^1/3 



(3.50) 



we are neglecting higher powers of Bjorken-x, which seems to be justified in the small-x limit. ^ 
However, it is easy to see from Eqs. ( |3.49| ) that if we expand the exponentials in them back into 
series, the series would be in the powers of 



Ql{x,A) AMV3 



We would then have 

W''{q+,q-,q 



0)~E 



Q' 



AMV3 
Q^x 



Q^x 



[l + dixA'/-' + dl{xA'/Y + ...] 



(3.51) 



(3.52) 



with c„'s and d^'s some x- and Q^-mdependent constants. We now have a problem: if we want to 
take a small-x/fixed-Q^ limit (which is the same as taking p"*" large in Eqs. (|2.4|) and ( p.5D ), the 
eikonal formulas (|3.49|) would receive order- 1 corrections. Namely, if we do power counting in x in 
Eq. ( p. 521) , we see that subleading order-x non-eikonal correction to the (n + l)st term is of the 
same order as the nth term in the eikonal series 



Cn+l 



AMV3 
Q^x 



n+l 



dl^^xA^/^ ~ c„ 



A^A^ 
Q^x 



(3.53) 



as each of them is of the order x ". The parametric equality can be simplified to 

A2 ^2/3 



Q' 



1. 



(3.54) 



"'Note that q^ a ^ a/lcoh, where Icoh is the coherence length of the projectile (particles the _R-current decays 
into) in the x^-direction; smallness of q~^ a means that the coherence length is much larger than the size of the 
proton/nucleus. 



We see that the non-eikonal corrections become important at Q 
Q2 defined by 



Q2 with the momentum scale 



(3.55) 



This scale is similar to the saturation scale identified in scattering a dipole on a shock wave in [67] 
(see Eq. (4.28) there) and also in [78]. The only difference between Q2 and the saturation scale 
found in [67] is a factor of V^ contained in the latter which appears not to be present in Q2 (A 
is the 't Hooft coupling constant). This factor of vA is inherently present in any calculation of 
the dynamics of a fundamental superstring as it is a prefactor in the Nambu-Goto action (see 
e.g. [76]). At the same time any AdS/CFT-based calculation of i?-current correlators in the large- A 
limit appears to be A-independent. To date we have not found a satisfactory explanation of this 
difference by a/A and suspect that it may be related to some more fundamental questions concerning 
AdS/CFT correspondence. We leave this question open for further research.^ 

A priori the scale Q2 from Eq. ( p. 55 ) is not the only scale one can construct out of non-eikonal 
corrections in Eq. (|3.52|) . Equating the kth non-eikonal correction to the {n + m)th term in the 
series of Eq. ( p.52| ) to the nth eikonal term yields^ 



Cji+m 






n+m 



d: 



n+m 



{xA'/')' 



~ Cn 



A^A'/' 
Q^x 



(3.56) 



giving a possible new scale at which non-eikonal corrections should become important: 



Q 



^2^1/3 



m,k 



X 



(a;Ai/3) 



l/3xfc/m 



(3.57) 



Adjusting positive integers k and m in Eq. ( |3.57|) one can get the scale Qm,k as (parametrically) 
close to Qi as desirable. Still, as fc,m > we have Qm,k < Qi- In fact Qm,k^ are a multitude of 
scales below Qi and both above and below Q2- Since Qm,kS are the scales at which non-eikonal 
corrections are important, we conclude that, at least with this a priori analysis, one can not trust 
the eikonal formulas ( p.49|) for Q ~ Qm,k < Qi- 

Therefore, the conclusion of our power-counting analysis is that the eikonal expressions in 
Eqs. ( |3.49|) are valid only at 



Q' >Ql{x,A) 



(3.58) 



^Indeed /i may depend on A, as was suggested in Appendix A of [82]: however this would not explain the difference 



between Eq. ( |3.55| ) and Eq. (4.28) in [67], as A-dependence in /x would modify both of them in the same way. 

^It may happen that one of these terms is real, while the other one is purely imaginary: if one then insists on 
equating only real terms to real terms and imaginary terms to imaginary terms, one should equate the nth eikonal 
term to the (2 A:)th correction in the {n + 2 m)th term, obtaining the same result as below. Powers of i in Eqs. ( 3.39| ) 
and (3.49) would insure that the terms we compare are either both real or both imaginary. 
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or, strictly speaking, only for Q^ ^ Q\{x,A). To clarify whether this really is the region of validity 
of the eikonal approximation or whether it may actually be applicable at Q^ < QKx, A) one has to 
find the exact expression for the hadronic tensor. 
Note that since 



-| K_ j^ 



Qm,k ~ Ql "" QS 



(3.59) 



the problem of R-current DIS on a shock wave remains a problem with only two momentum scales 
Qi and Q2. The series (p.52 ) can be written in terms of these two scales as 



n'^^(g+,g-,g = 0) ~ J] 



01 



2n " 



1 + d. 



I ( Q'i\ , j2 (Q^ 



" VQi 



+ dt 



Qi 



+ 



(3.60) 



The important conclusion of this Subsection is that the DIS process at strong coupling is a two-scale 
problem. 

3.3 Beyond the Eikonal Approximation: Perturbative Solution 

Let us construct the hadronic structure tensor perturbatively by exactly calculating the terms in 
the series of Eqs. ( ^.21 ) and ( 3.33| ) order-by-order in /i. We denote the nth term in each of those 
series by 



U^\q\q-,q = 0) 



A^S± Nl /i 



a »7r^ 



^'5.,Q'(g~)' I dx- I dy-e^^^^'^'-y-^ f dzz^KiiQ. 



X 



/^ 



-oo 

n— 1 



iq'Y G,{-i q- 9_) z^ e{x-) e{a - X-) z Ki{Q z) bix' - y~) 



(3.61) 



and 



r(") 



nL"l(g+,g-,g = 0) 



A^Si_ Nl li 



a 32 7r2 



-^Q^ I dx- I dy-e^'^^^'-y^ f dzz^KoizQ) 



/i 



la 
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with n = 1,2 Below we estimate the n = 1 and n = 2 terms. 



(3.62) 



3.3.1 Leading Order 

For the n = 1 term a quick calculation readily gives the leading-order (LO) terms 
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The same expressions would be obtained if one expands the eikonal formulas ( p.42| ) and (|3.43| ) to 
order-/i. Indeed the eikonal approximation of the previous Subsection only modifies the gauge field 
propagators sandwiched between the graviton exchanges in Fig. |1| along with the interaction vertices 
of Fig. 1^, since we only modify the Green function as shown in Eq. (|3.36|) and the longitudinal 
integrals over positions of graviton-gauge field vertices, as follows from Eq. ( |3.4CI| ). For n = 1 we 
only have one graviton exchange: the propagators to the left and to the right of the graviton-gauge 
field vertex (each of which giving Ki{Q z) in Eq. ( p.21| ) and Kq{Q z) in Eq. ( p.33|) ) are exact, even 
in the eikonal approximation, and the longitudinal integral ( |3.39|) is carried out exactly in this case. 
Therefore the exact one-graviton exchange results in Eq. ( p.63|) are the same as the order-/! terms 
in the eikonal formulas ( p.42| ) and ( |3.43|) . 

As we did above we replace fj, and a using 
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a 2x AV3 
and 
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2x 
Eq. ( 3. 63 ) can then be written in terms of x and Q^ as 
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n!]' = % A^ 5x ^ ^, n« = A^ S, -|L _^ ^. (3,66) 

■' 40 TT"^ x'' 240 TT^ [q y x'^ 

3.3.2 Next-to-Leading Order 

We start by analyzing the transverse components of Tl^^ . To find the transverse hadronic tensor at 
the next-to-leading order (NLO) we put n = 2 in Eq. ( p.6ip and employ the definition of the Green 
function in Eqs. (|3.9|) and (^.6]) obtaining 



a a oo 

nlf = 5^^^^, {^y Q'iQ-r jdx- I dy- e^'^' ^^-y-^ J dzz'MzQ) 



oo 

X f dz'h (z< ^-i2q-d^\ Ki {zy ^/-i2q-dS] z'^ Ki{z' Q) 5{x' - y~) (3.67) 



where we have used 5{x~ — y~) to replace 6{x~) 9{a — x~) by 0{y~) 9{a — y~) and employed the 
latter to modify the limits of ^"-integration. As before 9_ = d/dx~ and z^ (<) = max (min) {z, z'}. 
We write 

oo 

S^x--y-)= /^e-^('''+"H^"-^"). (3.68) 
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The +i e regulator is inserted to impose causality: it makes sure that 

1 
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.S(x -y ) = e{x -y ) 



(3.69) 



for 1/d- defined in Eq. ( p.8|) . Using Eq. (|3.68| ) we rewrite Eq. (|3.67|) as 
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Performing x~ and ?/~ integrations in Eq. ( p. 701) yields 
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We have defined 
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In arriving at Eq. ( p.71| ) we have also used the fact that Q < Q^ = — 2g+ g^, such that, since a; > 0, 
then, due to Eq. (|2.10|) , g~ > and we have g"*" < 0. 

The ^-integral in Eq. ( |3.71D is analyzed in Appendix ^, where the z and z' integrals are also 
carried out. The result is (see Eq. ( |A7| )) 
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If one performs ^/-integration in Eq. ( p.73| ) one obtains an answer expressed in terms of special 
functions. However this does not appear to make the expression ( |3.73| ) more transparent: we will 
leave it in the integral form. Eq. ( p. 731 ) is our exact result for the transverse components of the 
hadronic tensor 11 jj at the order /i^. 



To explicitly find corrections to the eikonal expression we expand Eq. ( p.73|) in powers of a and 
integrate over y to obtain 
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Using Eqs. ( p.50| ), ( p.64| ), and ( p.65| ), we rewrite Eq. ( |3.74| ) in terms of x and Q 
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Finally, employing Eqs. ( p.48| ) and ( p.55| ) we rewrite our result ( p.75| ) as 
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We now move on to the longitudinal components of the hadronic tensor. Eq. ( |3.62| ) gives 
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The rest of evaluation proceeds along the same lines as for the transverse components of 11^''. 
Similar to Eq. ( p.71| ) we write 
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For evaluation of z-, z'- and ^-integrals in Eq. (|3.78|) the reader is referred to Appendix ^. Using 
Eq. ( |A12| ) there we write 
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This is our final exact result for 11 at the order-fi^. 

Expanding Eq. ( p.79|) in the powers of a yields 
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or, in terms of Bjorken x and Q^ 
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and in terms of Qi and Q2, 
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(3.82) 



One can see that the form of the hadronic tensor suggested in Eq. ( |3.6CI| ) is exphcitly confirmed 
by our results in Eqs. ( p. 761 ) and ( |3.82| )! The prefactors of the square brackets of Eqs. ( |3.7(j| ) and 
( p.82| ) can also be obtained from the eikonal expression (|3.49|) . For Q"^ = Ql{A) the second terms in 
the square brackets of Eqs. (|3.76| ) and (|3.82|) become parametrically comparable to (and numerically 
much larger than) the leading-order results given in Eq. (p. 66 ). As we noted above, this indicates 
the breakdown of the eikonal formula (|3.49|) at Q^ ~ Q2{A). 



3.4 Brief Summary of Our Results and Expressions for Structure Functions 

Let us briefly summarize the results of this Section. We have written down exact general expres- 
sions ( p.21| ) and ( p.33| ) for the two independent components of the hadronic tensor 11^^. These 
expressions do not appear to be easy to evaluate in general since they involve multiple iterations of 
the Green function operators Gz and G^- Instead we have employed two approximations aimed at 
understanding the structure of the full solution. 

We flrst re-derived the components of the hadronic tensor in the eikonal approximation [65,69- 
71] obtaining 
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(3.83b) 
We can re-write this result in terms of dimensionless structure functions Fi and F2 deflned by 

^i(^'Q') = 7^. Imni(x,Q2), F2(x,g2) = ^^ Imn2(x,Q2) (3.84) 
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where we replaced the conventional proton mass by the typical transverse momentum in the shock 
wave A. For p = q = case considered here we have 



nl(g+, q-,q = 0) = 6'^ Uu n++(g+, q~ , q = 0) 
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Combining Eqs. (pD, j^M) and (^) we write [65,69-71]^ 
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(3.86b) 



where, as usual, F^ = F2 — 2 x Fi. 

We have argued that the eikonal expressions ( p.86| ) apply only for Q^ > QKxjA). To demon- 
strate this explicitly we have evaluated the hadronic tensor at the orders /i and /i^ going beyond 
the eikonal approximation. Our above calculations can be summarized as follows: 
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(3.87a) 



(3.87b) 



Clearly for Q^ ^ Qil^) ^^e third term in the brackets in each of the equations (|3.87|) becomes 
comparable to the first (leading-order) term in the series, thus generating an order-one correction 
to Eqs. ( p.83| ). Corrections to the structure functions F2 and F^ result from the imaginary parts 
of Hi and 112 in E<1- ( |3.87| ): at the order of the calculation shown in Eq. ( |3.87| ) such corrections 
appear to stem only from the last terms in the square brackets, which are always smaller than the 
2nd terms contributing to the structure functions. However, it is clear that an order-yU^ calculation 
would generate imaginary terms oc iQ\Q\/Q'^ in the square brackets of Hi and 112, which would 
become comparable to leading large-Q^ contributions to the structure functions for Q"^ ^ Q^iA), 
generating corrections to ( |3.86|) . Such order- /i^ calculation, while conceptually straightforward, is 
technically rather involved. We have verified that the terms oc ^QiQi/Q^ "^o indeed arise in such 
calculation: the determination of the exact numerical prefactors in front of such terms does not 



^Note that iVc-counting here agrees with the perturbative calculations of the DIS structure functions for partons 
in color-adjoint representation with the target nucleus made of nucleons with N^ valence partons each. 
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seem to be important for the conceptual conclusion about the breakdown of the eikonal formulas 
( p.86| ) at Q"^ ^ Qli^) which we draw here. 



However, the applicability region of Eqs. ( p.86|) is not simply Q'^ ^ Q^i^A). In fact, as was 



detailed in Sec. 3.2.2| , non-eikonal corrections to the higher-order terms in the eikonal series lead to 



the apphcability region of the eikonal approach being reduced to Q^ > Qi{x, A), as such corrections 
become important at scales arbitrary close to (but smaller than) QKxjA). At the moment we can 
not asses the net size and the effect of such corrections: this may require knowing the exact solution 
of the problem (i.e., the exact solution of Eqs. ( |2.14| )). 

4. Discussion of Momentum Scales in DIS 

Above we have shown that i?-current DIS on a shock wave of finite longitudinal extent at strong 't 
Hooft coupling is described by two momentum stales, Qi and Q2- This seems to be natural since 
the finite-size shock wave is described by two dimensionful scales: /i and a. 

Our conclusion also appears to be in qualitative agreement with the calculation performed 
in [67]: there DIS process on a shock wave was modeled by a quark-anti-quark dipole scattering on 
a shock wave. Indeed in QCD this is how DIS process takes place: virtual photon splits into a quark- 
anti-quark pair which then scatters on a target proton or nucleus (see e.g. [20, 21] and references 
therein). In [67] the dipole-shock wave scattering was described by calculating an expectation 
value of a fundamental Wilson loop in the shock wave background. For a shock wave of the type 
( p.l| ) which does not have any transverse coordinate dependence, the resulting forward dipole-target 
scattering amplitude N{r, s) was a function of the transverse dipole size r and the center-of-mass 
energy s. In order to obtain the DIS structure functions one has to convolute N{r, s) with the 
light-cone wave function of a virtual photon splitting into qq pair, which results in r being dual to 
1/Q (see e.g. [68,93]). The dipole amplitude N{r, s) from [67] had two momentum scales associated 
with it. One scale was defined by the condition (fi/a) r'^ = const and indicated the qq separation r 
at which the classical string solution became complex-valued. To translate this scale into Q^ and 
X variables we note that the calculation in [67] was done in the rest frame of the dipole. Taking 
into account Lorentz properties of fi/a we see that to generalize this condition we should replace r^ 
by {q~Y/Q^ (as the appropriately transforming projectile-related parameters), such that this scale, 
which we label Qs, is defined by the condition 
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The second scale describing DIS obtained in [67] was the scale Q2{A). In complete analogy with 
perturbative QCD calculations, defining saturation scale by N{l/Qs, s) = 1/2 (half of the black 
disk limit of A^ = 1) [41], the calculation in [67] obtained Qs ~ Q2{A.) (see also [78]). It was also 
observed that the dipole amplitude N{r, s) became independent of energy /Bj or ken x at high energy 
for a broad range of values of r both inside and outside of the saturation region. 

We see that the two scales Qi and Q2 we have obtained above were also present in the calculation 
of [67]. In the small-x regime, when x A^'^ <^ 1, the lower scale found in [67] was Q2 which 
corresponded to the saturation scale. The larger scale Q3 from Eq. ( ^4.2| ) can also be understood: 
in [67] the string solution considered was static, which is strictly speaking only valid for a shock wave 
of infinite extent. Hadronic tensor for such shock wave can be obtained from our exact equations 
( ^.21| ) and ( |3.33D by taking a — )■ 00 limit while keeping ft/a fixed: this would simply remove theta- 
functions, making the series in Eqs. ( p. 211) and ( p.33| ) a power series in (/i/a) {q~Y /Q^ ■ The problem 



now is described by only one scale — the scale Q3 from Eq. ( [1.2|) . 

To reconcile this single-scale result with our two-scale conclusion above one could argue follow- 
ing [65] that for DIS on an infinite-extent shock wave what matters is the size of the interaction 
region between the /^-currents and the shock wave. One should therefore replace the shock wave 
longitudinal width a ~ A^^^ /p^ by the typical longitudinal separation between points x~ and y~ 
in the i?-current correlator (|2.6| ): the latter is ~ l/l'?"''! = l/(xj9"'"). Hence one has to replace 

^1/3 ^ i. (4.3) 

X 

Under such replacement both scales Qi and Q2 become equal to Q^ and the problem becomes 
single-scale. This is also why the saturation scale Qs = Q2 found in [67] is in agreement with the 
results of [65,66] for the infinite medium: under the substitution (|4.3|) one has (^2)^ = A^A^'^ — >■ 
A^/x^ = (Qs)^. (The discrepancy by a factor of V^ with A the 't Hooft coupling that we mentioned 
above still remains indeed: we can not explain it at the moment.) 

Let us now return to the finite-extent shock waves. While both our above analysis and the 
calculations of [67] have the same conclusion about DIS at strong coupling being a two-scale problem, 
one may still worry about the physical interpretation of the scales Qi and Q2 we found. The 
calculation in [67] considered DIS on a finite-size shock wave but, as a first approximation, employed 
the static limit of the string configuration, strictly speaking valid for an infinite shock wave only: the 
question whether in scattering on a large but finite shock wave the string has enough time to quickly 
settle onto its static configuration remains to be answered (see [70] for the analysis of the problem 
for a thin shock wave). Moreover, building on the analogy with the complex trajectory method 
in Quantum Mechanics the classical string solutions found in [67] were analytically continued into 
the complex-valued domain. Justification of such a procedure may be needed in the string theory 
context. Therefore we will try to discuss the possible physical meanings of Qi and Q2 with an open 
mind while temporarily ignoring the pre-existing results of [67]. 
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Indeed the physical meaning of the scales Qi and Q2 would have probably been more manifest 
if the exact result for F2 structure function was known. Instead we have the eikonal expression 
( pGal ) due to [65,69-71], which is valid for Q"^ > Ql{x, A). At large Q^ > Ql{x, A) it gives [69-71] 



i^2(x,Q' 



while at small Q^ it reduces to [69,71] 



2»Q2(x,A) 



F2{x,Q' 



\Q'^<€.Qfix,A) 



5, 



^1 



i8iy|AMV3 

35 7r3 x2Q2 ■ 



(4.4) 



N^ Q2 g2(^^^) 



64 7r3 Ai/3 



In 



Q' 



(4.5) 



though it is not clear how reliable Eq. (^4.5|) is in light of Q^ > Qi{x,A) applicability constraint of 
Eq. (|3.86a|) . We will proceed under assumption that Eqs. (|^ and ( [4. 51) are quahtatively correct, 
i.e., that F2{x,Q^) has a maximum at Q^ ~ Ql{x,A) and it decreases as Q^ becomes either larger 
or smaller than QKx, A). 

The exact effect of the scale Q2{A) on F2 is not clear from the above calculations and it appears 
that to clarify it one needs to solve the problem exactly. In the meantime we argue that Eqs. ( p.87|) 
indicate that the structure functions would change quite significantly at Q = Q2{A). We can only 
guess the exact effect of Q2 on F2. If we believe that F2 already decreases with decreasing Q^ for 
Q < Qi, as seems to follow from Eq. ( [4. 5] ) which we choose to believe at least at the qualitative 
level, and combine this with the fact that, on general grounds, F2 should go to for Q^ -^ 0, we 
conclude that it is probable that for Q < Q2 the structure function F2 would continue to decrease 
with decreasing Q^, probably decreasing faster than it was for Q2 < Q < Qi- The sketch of our 
guess/tentative understanding of F2(x, Q"^) is shown in Fig. ||. 




Q. 



q; 



Q^ 



Figure 2: A sketch of the F2 structure function of i?-current DIS at strong coupling as a function of Q^ 
based on our understanding/guess of the exact AdS/CFT prediction (see text). 
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Assuming that the sketch in Fig. ^ accurately represents the structure function F2 given by the 
exact solution of the i?-current scattering problem in AdS/CFT, we propose the following three 
possible interpretations of the physical meaning of the scales Qi and Q2- 



1. The first, and, in our view, the most probable option, stems from comparing F2 structure 
function in Fig. |^ to what one has in QCD at small coupling and/or to the actual data reported 
by DIS experiments. The prediction of CGC/saturation physics is that the F2 structure function 
scales as [21,94-96] (for a review see [27])® 

Ff^^ (xQ\ Q^<^Ql (4.6a) 

Ff «^ oc {QT'', Q' > Ql (4.6b) 

It is important to note that in the small- a; CGC/saturation physics framework F2 structure function 
never decreases with increasing Q^\ Therefore it may seem hard to reconcile the decrease of F2 with 
Q^ shown in Eq. ( |4.4| ) with small- a; CGC/saturation physics. Moreover, if one remembers that in 
weakly-coupled QCD F2 is given by the sum of quark distributions xq{x, Q"^) over all flavors the 
decrease with Q^ may seem even more puzzling: if we use the standard (albeit somewhat simplified) 
interpretation of xq{x, Q"^) as the number of quarks at a given value of Bjorken x with transverse 
momenta kx < Q, then it would appear that xqi^x^Q"^) along with ^2(0;, Q^) can never decrease 
with Q^, since the number of quarks with kx <Q can only increase with Q^ . 

However such arguments are not entirely correct. It relies on a simple perturbative relation 
between F2{x,Q'^) and distribution functions, which may be modified at strong coupling. On top 
of that^, at large-Q^, when the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
[98-100] is important, the relation between, say, the gluon distribution function xG{x, Q^) and the 
unintegrated gluon distribution is not simply 



xG{x, Q'') = / dkl 0(x, 4) (4.7) 

but, instead, is given by [101-105] ^^ 



xG{x,Q') = / dk^(l){x,k^,Q'). (4. 



The power of Q^ in Eq. (4.6b|) is given by the fixed-couphng approximation. This power changes when running 



couphng corrections are included [41,97]. For our purposes we only need the power to be positive and smaller than 
1, which is true at both the fixed and running coupling. 

^We would like to thank Genya Levin for pointing out this argument to us. 

^°0f course Eq. (4.8) should get substantially modified and ceases to be valid also at low-Q^ inside the saturation 



region {Q < Qs) due to multiple-rescatterings, non-linear evolution, and other higher-twist corrections. 
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As usual for distribution functions, Q^ is the renormalization scale: in the spirit of the leading 
logarithmic approximation we put it as the upper cutoff on the k'^ integral in Eq. ( [4.8|) . Eq. ( |4.8| ) 



shows that when one goes beyond the leading logarithmic small-x evolution approximation, and 
includes DGLAP evolution [98-100] as well, the unintegrated gluon distribution 0(x, /cf., Q'^) itself 
becomes a function of Q^. It implies that the highest transverse momentum of a "real" parton in 
the proton's wave function is kx-, while the wave function is evolved using DGLAP to the scale Q^, 
such that the evolution from k"^ to Q^ is due to virtual corrections to DGLAP only, resulting in a 
form-factor in the definition of 0(a:, /c|., Q'^) (see [102-105] for more details and for the definition of 
0(x, /cf., Q^)). Now, as Q^ increases, the unintegrated distribution function (l){x,k^,Q'^) decreases, 
as the probability of no real gluon emissions between k^ and Q^ decreases with increasing Q^. It 
is thus possible that at large enough Q^ this decrease with Q^ in would dominate in Eq. ( |4.8|) , 
resulting in the gluon distribution function xG{x, Q"^) decreasing with Q^. (The same argument can 
be applied to quark distributions.) 

It is important to point out that if one wants to interpret (j){x, k^, Q^) in Eq. ([4.8|) as the number 
of gluons with transverse momentum kr, this number would depend on the momentum of the probe 
(or, equivalently, on the renormalization scale) Q. The reason behind this Q-dependence is that 
0(x, k^, Q'^) really gives the number of gluons at kx with the condition that there are no gluons with 
higher transverse momenta in the hadronic wave function. It appears to be impossible in general 
to define unintegrated gluon distribution independent of Q^, which would simply give the number 
of gluons at fcy without any exclusive conditions. Thus the probabilistic interpretation of the gluon 
distribution xG{x, Q'^) as the number of gluons with fcy < Q is not valid once the full DGLAP 
evolution is included. (Again the same applies to quark distributions.) This is why the falloff of F2 
with Q^ presents no contradiction. 

To visualize how a distribution function (and therefore a structure function) may decrease with 
Q^ and to determine at what Q^ these functions start decreasing let us consider a simple but 
realistic toy model. ^^ Take the gluon distribution given by the solution of the leading- logarithmic 
fixed-coupling DGLAP evolution equation: 



b+i 00 



xG{x,Q'')= J ^x"(^l G^iQl). (4.9) 

b—i 00 

Here b is an arbitrary real number and Qo is the initial scale of DGLAP evolution. For simplicity 
we assume that there are no quarks in the toy theory we consider. We also assume a particularly 
simple toy form of the gluon-gluon splitting function [106] 

iGGico) = ^^(--1). (4.10) 

n \UJ J 



^We are grateful to Genya Levin for this argument as well. 
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This splitting function has the correct residue of the small-x pole at w = 0. The term (—1) in the 
parenthesis of Eq. ( [4.1(]| ) mimics all the non-small- a: terms in the actual splitting function. It also 
makes sure that the momentum sum rule 

7gg(1) = (4.11) 

is satisfied. 

At small-x and large-Q^ the integral in Eq. ( |4.9| ) can be evaluated in the saddle-point approxi- 
mation with the saddle point at 






and the gluon distribution given approximately by 



Q' 



2\ 7Gg('^sp) 



q: 



xG{x,Q') ~ x"^^'' {-^] . (4.13) 



This distribution function is a decreasing function of Q^ for •yGci^sp) < 0, which means Usp > 1. 
Therefore the gluon distribution decreases with Q^ for 



(1 \ as Nc 
x) ■ (^-^^^ 

This is indeed a very large scale for small- a; , but for larger-x it becomes small enough for decrease 
of F2 with Q"^ to be seen experimentally at HERA. Note that at large 't Hooft coupling the scale 
in Eq. ( |4.14 ) is not necessarily large. 



Eq. ( [4. 141) illustrates a known fact that at very large Q^ distribution functions (and, therefore, 
structure functions) do decrease with Q^ even in the perturbative picture. Therefore, combining 
this result with Eqs. ( |4.6|) we now see that in perturbative QCD the F2 structure function looks 
qualitatively as shown in Fig. ^ if we identify the scale Q2 with the saturation scale Qg and the 
scale Qi with the scale Q^^cr i^om Eq. ([4.14|) at which F2 starts falling off with Q"^. 

Therefore our first guess at the physical meaning of Q2 and Qi is to identify them with Qg 
and Qdecr correspondingly. In [65,66,71] it is shown that the scale Qi is essential for satisfying the 
momentum sum rule: this seems to confirm our conclusion since Qdecr results from satisfying the 
same momentum sum rule of Eq. (|4.11| ). It is possible that at strong 't Hooft coupling, just like 



at small coupling ««, energy conservation effects come in at a different Q^-scale from unitarization 
effects. 

2. The second interpretation of Q2 and Qi we propose is to leave the interpretation of Q2 as the 
saturation scale, but to suggest that Qi is the extended geometric scaling scale kgeom [95,96,107]. 
Extended geometric scale kgeom is the scale such that for Q < kgeom the structure functions are 
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functions of Q/Qs only [95,96]. In CGC usually kgeom > Qs [96], which supports our hypothesis 
here. Also, if we accept Eq. ( [4.5|) as being at least qualitatively correct for the exact AdS/CFT 
prediction for F2, we can see that F2 is almost completely x-independent below Qi, and probably 
can be written as a function of Q/Q2, which also supports the suggestion that Qi could be kgeom, 
since in CGC F2 is a function of Q /Qs for Q < kgeom [96]. Indeed the relation between kgeom and Qs 
has to be modified at strong coupling in comparison to the weak-coupling CGC result [27,96,108]. 

The main problem with this scenario is that in CGC for Q > kgeom the structure function F2 
keeps increasing with Q [27,96], while in AdS/CFT calculations one obtains F2 decreasing with Q 
for Q > Qi as one can see from Eq. (|4.4| ) and from Fig. 0. Therefore our second hypothesis does not 
seem to be in agreement with the shape of the plot in Fig. |^ for Q > Qi, which makes it somewhat 
less compelling than the first one. 

3. Finally one may accept the viewpoint advocated in [65,66,70,71] and identify Qi with the 
saturation scale. Indeed the similarity between Eqs. ( [4.6a] ) and ( ^75| ) seems to suggest that this is 
correct. However, as we argued above, Eq. (|4.5|), while obtained by eikonal methods, lies outside 
the region of applicability of the eikonal approximation and should be questioned. Also in CGC 
the structure function F2 continues growing with Q^ for Q > Qs, as one can see from Eq. ( [4.6b|) , 
in disagreement with Eq. (|4.4|) , casting more doubt on this third possible scenario. One should 
also mention that x-independence of structure functions was observed at large coupling in [67] for 
Q > Qs'- therefore x-independence of Eq. (|4.5D may not yet signal saturation. 



An important question remains regarding the physical role of the scale Q2. In traditional CGC 
literature there are no important scales below Qs- One may speculate that Q2 may be the scale 
at which other higher twist effects, such as pomeron loops, may become important (see e.g. [20]). 
While possible in principle we believe further research is needed to test this assumption. Pomeron 
loops are suppressed by powers of A, while the scale Q2 does not have any A- suppression compared 
to Qi, exhibiting the opposite A-enhancement. In principle, until the exact solution of the problem 
is found, it may also be possible that nothing of physical importance happens at the scale Q2, 
though such conclusion is hardly likely, since a whole class of terms becomes important at this 
scale, as one can see from Eq. ( |3.60| ). The scale Q2 is known to play an important role in heavy 
ion collisions modeled in AdS/CFT: as was shown in [79, 82] in a strongly-coupled collision shock 
waves stop at the light-cone time xj^p ~ ^/Q2 in the center-of-mass frame. It is probable that the 
scale that determines the stopping time in a shock wave collision should play some role in DIS as well. 



Indeed an exact solution of the i?-current DIS problem is needed to conclude whether one 
(if any) of the above-listed possibilities is correct. Unfortunately an exact analytic evaluation of 
Eqs. (|3.21|) and (p.33|) appears to be a rather difficult problem at present. 
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A. Some useful integrals 

We start by integrating over ^ in Eq. (|3.71|) , namely we need to find 



Ro = ^ 



d^ 1 

2n (1-e + ^e) 



2 _ g'ii3+a(l-C+*<:) _ g-ig+a(l-^+Je) 



X 



oo oo 

I dzz'^K^izQ) I dz' h (^< Q Ve-^e) Ki {z> Q V^-^e) z"' K^ {z' Q). ( Al) 



Using the series representations of the modified Bessel functions 



^i(-) = E 



m=0 

oo 



K,{z) = J2 



rfm)rfm + l) V2 



rra=0 



r(m)r(m + l) V2 



2m-l 



2m~l 



In ( ^ ) - 7^ H^) 

2/ 2m 



(A2a) 
(A2b) 



we see that /i (2< Q y^^ ~ ^ ^ ) -^i (^> Q V^ ~ ^ ^) from Eq. ( |A1|) has a branch cut discontinuity for 
^ G (— oo + ie,0 + ie]. The complex structure of the integrand in Eq. (|A1|) is depicted in Fig. ^ 
The integrand has a branch cut we have just mentioned, along with a possible pole at ^ = 1 + i e. 
While strictly speaking there is no pole at ^ = 1 + i e in the full expression in Eq. (|Al|), individual 
terms in the square brackets in Eq. ( |A1| ) lead to contributions to the integrand containing this pole. 
Since q~^ < 0, the last term in the square brackets of Eq. ( |A1[ ) demands that the ^-integration 
contour be closed in the lower half-plane: since there are no poles or branch cuts in the lower 
half-plane, we can discard this term. For the first two terms in the square brackets of Eq. ( |Al| ) we 
have to close the integration contour in the upper half-plane (for the very first term the direction 
of contour closing is actually dictated by the large-argument asymptotics of the modified Bessel 
functions). Picking up the pole at ^ = 1 + ie and wrapping the contour around the branch cut 
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Im^ 



Re^ 



Figure 3: The complex ^-plane structure of the integrands of Eqs. (Al) and ( |A8[ ) (see text) 



yields 



R-2 



d^ 



2 J (1-0^ 



2 _ g«9^a(l-?) 



-oo 
oo 



dzz''Ki{zQ)Ji UQV-^ 



X 



f dz z"- K,{z Q) I dz' h (^< Q Vt) Ki (2> Q v^) z'' K,{z' Q) 





where we have used 



d 



C=i 



2 _ g«9+a(i-5) 






which can be inferred from Eqs. ([A2| ) . 

The first term on the right-hand-side of Eq. ([A3|) is straightforwardly evaluated as 



di 



2 J (l-eP 

— oo 





2 — g*9"'''^(l-?) 
-4j- L 

i 192^ /• d^ 



dzz^Ki{Qz)Ji (zQ^/^ 



2 Q12 J (1-0' 

— oo 
oo 

i 1922 r dy 



2 — g^e+aCi-^) 



^(1 + 



2 gi2 y (i + y)i2 





y(i-z/)' 



(1-0^ 

2 _ g^'J^^Ci+s/) 



(A3) 



(A4) 



(A5) 



where y = — ^. The second term on the right-hand-side of Eq. ( |A^) is a linear polynomial in a, 

# + #'a, (A6) 



with the coefficients depending on Q. We know from the eikonal approximation (see Eq. (|3.40| )) 
that i?2, if expanded in a series in the powers of a, should start at the order a^. The same can be 
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inferred from Eq. (Kl), though we note that a power-series in a expansion in the integrand there 
gives finite results only at the order a^, not allowing to learn anything about higher powers of a. 

Requiring that the series in powers of a for R2 starts from a^ along with Eq. ( |A6| ) shows that 
the second term on the right- hand- side of Eq. (X3) simply cancels the constant and linear in a terms 
in the first term on the right-hand-side of Eq. (|A^). (We have also checked by explicit numerical 
integration that this is true.) Adding extra terms to remove the constant and linear in a terms in 
Eq. ( [A5| ) we obtain our final answer for R2: 



R2 



i 192^ 



dy 



2 Q12 J (1+^)12 




y{l-yf [l + zg+a(l + 2/)-e^''''"(^+^) 



(A7) 



We now need to evaluate 



i?2 



di 



2 _ g«9"'"a(l-C+ie) _ g-i(/+a(l-5+Je) 



27r (l-^ + ie)2 

—00 

00 00 

X I dzz'' Koiz Q) I dz' Jo {z< Q Ve-^e) ^0 (^> Q V^-^e) ^" ^o(^' Q)- (A8) 


needed for calculation of the longitudinal components of the hadronic tensor in Eq. (|3.78| ). We 
begin by employing the series representation for modified Bessel functions Iq and Kq 



Iniz) 



Kn(z) 



771=0 

00 

-E 



2m 



[r(m + l)]2 V2, 

1 /2x2m 



In ( - j — 'ip(m + 1) 



(A9a) 
(A9b) 



- [r(m + l)]2 V2. 

to infer that the complex ^-plane structure of the integrand in Eq. (|A8|) is the same as shown in 
Fig. 1^ above. Picking up the pole at ,^ = 1 + z e and wrapping the contour around the branch cut 
yields, similar to Eq. 



Xtn 



d^ 



2 J (1-0^ 

-00 

00 



2 — e* 9+ 1(1-5) 



dzz''KoizQ)Jo{zQ^/-^ 



X 



I dz z'' Ko{z Q) j dz' h (^< Q Vt) Ko (2> Q Vc) z" Ko{z' Q) 



d 



?=i 



2 _ g^g"*" «(!-?) 



(AlO) 



where we have used 

= T,ie{-i).h[zQ^^) j^(^z'Q^^) (All) 
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which follows from Eqs. ( |A9|) . The second term on the right-hand-side in Eq. ( A10| ) is again a linear 
polynomial in a, with the coefficients that can be fixed by requiring that the Taylor expansion of 
R2 in powers of a starts from order-a^. Imposing this condition and integrating over z in the first 
term on the right-hand-side of Eq. ( |A1CI|) we arrive at the final result 



lin 



i 642 



2Q 



12 



dy 



a + y) 



12 



[l-4y + y 



2\2 



1 + iq'^ a{l + y) — e 



iq+ a{l+y) 



(A12) 



where, as before, y = — ^. 
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